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ABSTRACT Cj 

An analytic solution for high-energy heavy ion transport assuming 
straightahead and velocity conserving interactions with constant nuclear cross 
reactions is given in terms of a Green's function. The series solution for the Green's 
function is rapidly convergent for most practical applications. The Green’s function 
technique can be applied with equal success to laboratory beams as well as to galactic 
cosmic rays allowing laboratory validation of the resultant space shielding code. 
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Introduction 


Several approaches to the solution of high energy heavy ion propagation 

have been developed (1-12) over the last 20 years. All but one (5) have assumed the 

straightahead approximation and velocity conserving interactions. Only two 

(5, 8) have incorporated energy-dependent nuclear cross sections. The approach by 

Curtis et al. (1) for a primary ion beam represented the first generation secondary 

fragments as a quadrature over the collision density of the primary beam. Allkofer 

and Heinrich (2) used an energy multigroup method in which an energy- 

independent fragmentation transport approximation was applied within each 

energy group after which the energy group boundaries were moved according to 

[_dE.l 

continuous slowing down theory I . Chatterjee et al. (3) solved the energy- 
independent fragment transport equation with primary collision density as a source 
and neglecting higher-order fragmentation. The primary source term extended only 
to the primary ion range from the boundary. The energy-independent transport 
solution was modified to account for the finite range of the secondary fragment ions. 
Wilson (4) derived an expression for the ion transport problem to first order (first- 
collision term) and gave an analytic solution for the depth-dose relation. In ref. 5, 
the more common approximations used in solving the heavy ion transport 
problem were examined. The effects of conservation of velocity on fragmentation 
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and the straightahead approximation are found to be negligible for cosmic ray 
applications. Solution methods for representing the energy-dependent nuclear cross 
sections are developed. Letaw et al. (6) approximate the energy loss term and ion 
spectra by simple forms for which energy derivatives are evaluated explicitly (even 
if approximately). The resulting ordinary differential equations in position are 
solved analytically similar to the method of Allkofer and Heinrich (2). This 
approximation results in a decoupling of motion in space and change in energy. In 
Letaw's formalism, the energy shift is replaced by an effective attenuation factor. 
Wilson adds the next higher-order (second collision) term (7). This term was found 
to be very important in describing 670 MeV/amu 20 Ne beams. The three-term 
expansion of ref. 7 was modified to include the effects of energy variation of the 
nuclear cross sections (8). The integral form of the transport equation (5) was 
further used to derive a numerical marching procedure to solve the cosmic ray 
transport problem (9). This method can easily include the energy-dependent 
nuclear cross sections within the numerical procedure. Comparison of the 
numerical procedure (9) with an analytical solution to a simplified problem (10) 
validates the solution technique to about 1 percent accuracy. Several solution 
techniques and analytic methods have been developed for testing future numerical 
solutions to the transport equation (11, 12). More recently, an analytic solution for 
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the laboratory ion beam transport problem has been derived assuming a 
straightahead approximation, velocity conservation at the interaction site, and 
energy-independent nuclear cross sections (13). 

In the above overview of past developments, the applications split into two 
separate categories according to a single ion species with a single energy at the 
boundary versus a broad host of elemental types with a broad continuous energy 
spectrum. Techniques requiring a representation of the spectrum over an array of 
energy values require vast computer storage and computation speed for the 
laboratory beam problem to maintain sufficient energy resolution. On the other 
hand, analytic methods (5, 9) are probably best applied in a marching procedure (9) 
which again has within it a similar energy resolution problem. This is a serious 
limitation since we require a final HZE code for cosmic ray shielding which has been 
validated by laboratory experiments. In the present paper, we examine new 
methods which appear to overcome these difficulties. 

The HZE Transport Problem 

In moving through extended matter, heavy ions lose energy through 
interaction with atomic electrons along their trajectories. On occasion, they interact 
violently with nuclei of the matter producing ion fragments moving in the forward 
direction and low energy fragments of the struck target nucleus. The transport 
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equations for the short-range target fragments can be solved in closed form in terms 
of collision density (ref. 5). Hence, the projectile fragment transport is the 
interesting unsolved problem. In previous work, the projectile ion fragments were 
treated as if all went straightforward. We continue with this assumption herein, 
noting that an extension of the beam fragmentation model to three dimensions is 
being developed (ref. 14). 

With the straightahead approximation and the target secondary fragments 
neglected (ref. 5), the transport equation may be written as 


i_ ± 

dx 3E 


Sj (E) + a,- 


<frj (X, E ) = X m jk <*k 4»k (*. E ) 
k > j 


( 1 ) 


where (x, E) is the flux of ions of type j with atomic mass Aj at x moving along the 
x-axis at energy E in units of MeV/amu, Oj is the corresponding macroscopic nuclear 
absorption cross section, Sj (E) is the change in E per unit distance, and mj k is the 
fragmentation parameter for ion j produced in collision by ion k. The range of the 
ion is given as 


Sj (F) 

The solution to equation (1) is to be found subject to boundary specification at x = 0 



and arbitrary E as 
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(0, E) = Fj (E) (3) 

Usually Fj (E) is taken as an incident laboratory ion beam spectrum or the cosmic ray 
spectrum. 

It follows from Bethe's theory that 


A p Zj2 _ 

S : (E) = — H_3 — s (E) 

J AjZp 2 p 


(4) 


for which 

Z i 2 V (5) 

Rj (E) = Rp (E) 

The subscript p refers to proton. Equation (5) is quite accurate at high energy and 
only approximately true at low energy because of electron capture by the ion which 
effectively reduces its charge, higher-order Born corrections to Bethe's theory, and 
nuclear stopping at the lowest energies. Herein, the parameter Vj is defined as 

Vj-Z^Aj (6) 

so that 


Vj Rj (E) = v k R k (E) 


(7) 


Equations (6) and (7) are used in the subsequent development, and the energy 
variation in Vj is neglected. The inverse function of Rj (E) is defined as 

E = R; 1 [ Rj (E) ] (8) 

and plays a fundamental role subsequently. 


Impulse Response 

One form of Green's function is the impulse response corresponding to a 8- 
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like source term at the boundary. We therefore seek a solution of 


5x 


- — SjCEl + Oj 

9E J J 


G jM (x, E; E’) = X m jk <*k G kM ( x , E ; E’) 

k 


(9) 


where the boundary condition is 

G jM (0, E; E’) = 5 jM 5 (E - E') (10) 

for which transport solution may be written as 

<t>j (x, H) = X fo Gjm (x, E; E’) F M (E’) dE' U 1} 

M 


The solution to equations (9) and (10) is straightforward, even if tedious (4, 7, 13) and 
is arrived at using the method of characteristics (15). The solution is expressed as a 
series as 

G jM (x,E;E') = X G$(x,E;E’) (12) 

i 

where 



(x, E; E’) 


= s ^ exp( -° )X)5jM5 


x + Rj (E) - R m (E ) 


(13) 
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(x, E; E') = rr- 1 — m jM a M 
S;(E) JM 







x + Rj (E) + t)' 


} 


so long as 

M R M (E) - x ] < Rj (E) < ^ R„ (E) - x 


where 


TV = ^ m ~R M (E')-^^- 


V M -Vj V M -Vj[ 


Rj (E) + x 


( 14 ) 


( 15 ) 


( 16 ) 


Otherwise Gj^j (x, E; E') is zero. After a complicated but straightforward 
manipulation, a similar result may be obtained g|^ (x, E; E') as 



(x,E; E’) = X- 

k Sj 


g jk g kM v j 

(E) I V M - v k | 


e 


* g M X M1 _g k x kl ' g j x jl 

-e 


" g M x Mu ' g k x ku * g j x ju 


( 17 ) 


where xm u , x ku , xm,, xjq are values of x M and x k evaluated at the corresponding 


upper and lower limits of x . and 


X M = 
x k = 


v M R m (E) - v k (R k (E) + x) + (v k - Vj) Xj 
v m (Rm (E) + x ) - Vm Rm (E) - ( v m * v j) x j 


/( v M' v k) 

/(v M - v k) 


( 18 ) 

( 19 ) 
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The requirements that x M and x k be bounded by the interval 0 to x-Xj yields 


1 

0 


X 


VI 

yT 

VI 


1 v k [ (E) + x ] - V M R m (E ) 


V M [ R M (E) + X ] - v M R m (E ) 

1 

< 

rr 

i 


V M-Vj ) 


as the appropriate limiting values in equation (17) when v M > v k > Vj . In the above 
brace, we always choose the most restrictive value for the limit. The requirement of 
equation (20) also implies the result that 



( 21 ) 


as the range over which the result of equation (17) is not zero. In the event that 

v k >v M> V hen 


0 




V M [ R.Vt (E) + x ] * V M R.M (E) 


V M-Vj 


< Xj < 


( 22 ) 


v k [ R k (E) + x ] - v M R m (E ) 


v k - Vj 


As a result of equation (22) 


Rk Itt r m <E)-x}SESRil{ R m (E)-^;x 

v k 


V M 


( 23 ) 


In the event that v M > Vj > v k , it follows that 
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0 - X J " V M [ R M (E) + x - R m (E)]/(v m - Vj) 
, v m R M (E)-v k R k (E)-v k x]/ (Vj - v k ) 


where the lesser of the three values in the brace is used as the upper limit of x- for 
which G (2) of equation (17) is not zero. As a result of equation (24) 


Rm L R M (E) - x j < E < R k ^R m (E')-x 


Note that 


a . ( v k ' v j) _ ( v M* v j) 


Higher-order terms are similarly derived. Approximate expressions have been 
obtained as 

G$ (x, E; E) - I "jk<V.. a i"-2M ea.k.ji"-.j„. l M)/(E iu -E jl ) (27) 

tJl,- Jn-2 


-1 V M 


{r m (E) - x) 


R M (E) -X 


{rm (E)-X. 


V M >V k >Vj 
V k>V M >Vj 
v M > Vj > v k 


R k V M R M OB) * v j X /V k 


Rm Rm (E) • v j */ v m 


R k ' tt r m(e')-x 

I v k I 


V M >V k >Vj 

v k > v M > Vj 


V M > Vj > V k 
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and the g-functions of n + 1 arguments are defined as 

g 0"i) = exp ( - 0j t x) (30) 

„ r, \ ; ; n g (jl* J*2 • ' •> jn-1. jn) * g (jl> J2 ’ ' jn - 1 » jn + l ) (31) 

g 0l.j2.-”jn.jn + l) = 

The expression for given by equation (27) is taken as zero unless 

% < E < E ju (32) 

Portions of the Green's function are shown for incident 20 Ne beams at 
E' = 600 MeV/amu at x = 20 cm in figures 1 to 3. The contribution from G (1) is 
shown as the solid curve, G (1) + G (2) is shown as the dash-dot curve, and 
G (1) + G (2) + G (3) is shown as the dashed curve. The dash-double-dot curve 
representing the inclusion of G (4) in the sum can hardly be distinguished signifying 
convergence to a high degree of accuracy. A fuller presentation of the Green's 
function for 16 0 fragments is given in figure 4, and a presentation of the 12 C Green’s 
function is given in figure 5. From the present result, the solution for any arbitrary 
boundary condition may be found using equation (11). 

Although the present formalism presents a closed form solution for the more 
common form of the HZE propagation problem, many tasks remain before the HZE 
propagation problem is adequately solved. The inclusion of energy-dependent 
nuclear cross sections is known to be very important in obtaining accurate solutions 
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to some problems (16). Treating the momentum spread of the fragments is more 
complicated for the higher-order terms. The inclusion of the light fragment spectra 
is a difficult challenge (17). Finally, the three-dimensional aspects of the problem 
have only partially been treated (14). Even these shortcomings of the HZE 
propagation problem remain without the mention of uncertainties in nuclear cross 
sections (18) or atomic/ molecular cross sections. Clearly much work remains. 
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FIGURE CAPTIONS 


Fig. 1 - Sequence of approximations of the 17 0 flux spectrum after 20 cm of water 

for first order (-), second order (- • -), third order (—), and fourth order ( ) 

theories. 


Fig. 2 - Sequence of approximations of the 16 0 flux spectrum after 20 cm of water 

for first order (-), second order ( — ), third order (-), and fourth order ( ) 

theories. 

Fig. 3 - Sequence of approximations of the t2 C flux spectrum after 20 cm of water 

for first order (-), second order ( — ), third order ( — ), and fourth order ( ) 

theories. 

Fig. 4- Green’s function for 16 0 flux spectrum response to a 600 MeV/amu 20 Ne 
flux at the boundary. 

Fig. 5 - Green's function for 12 C flux spectrum response to a 600 MeV / amu 20 Ne 
flux at the boundary. 
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